Phylogenetics and Population Genetics of the Asian House Shrew, Suncus murinus-S. montanus Species Complex, Inferred From Whole-Genome and Mitochondrial DNA Sequences, with Special Reference to the Ryukyu Archipelago, Japan.

The house shrew (Suncus murinus-S. montanus species complex) colonized regions across southern Asia and the Indian Ocean following human activity. The house shrew is distributed on islands of the Ryukyu Archipelago, the southernmost part of Japan, but the evolutionary history of the shrew on those islands and possible associations between these populations and humans remain unknown. In this study, we conducted phylogenetic and population genetic analyses based on both nuclear and mitochondrial genome sequences of house shrews. Phylogenetic analyses based on mitochondrial cytochrome b (cytb) sequences revealed that shrews from the Ryukyu Archipelago showed strong genetic affinity to Vietnamese and southern Chinese shrews. Demographic analyses of cytb sequences indicated a rapid population expansion event affecting the haplotype group in Vietnam, southern China, and the Ryukyu Archipelago 3300-7900 years ago. Furthermore, gene flow between Ryukyu (Yonaguni Island) and Taiwan and between Ryukyu and Vietnam inferred from f4 statistics of the nuclear genomes suggested repeated immigration to Ryukyu in recent years. The present study demonstrates that the Nagasaki population has a different origin from the Ryukyu population. These findings elucidate the complex pattern of genetic admixture in house shrews and provide insights into their evolutionary history.

Inference of selective forces on house mouse genomes during secondary contact in East Asia.

The house mouse (Mus musculus), which is commensal to humans, has spread globally via human activities, leading to secondary contact between genetically divergent subspecies. This pattern of genetic admixture can provide insights into the selective forces at play in this well-studied model organism. Our analysis of 163 house mouse genomes, with a particular focus on East Asia, revealed substantial admixture between the subspecies castaneus and musculus, particularly in Japan and southern China. We revealed, despite the different level of autosomal admixture among regions, that all Y Chromosomes in the East Asian samples belonged to the musculus-type haplogroup, potentially explained by genomic conflict under sex-ratio distortion owing to varying copy numbers of ampliconic genes on sex chromosomes, Slx and Sly Our computer simulations, designed to replicate the observed scenario, show that the preferential fixation of musculus-type Y Chromosomes can be achieved with a slight increase in the male-to-female birth ratio. We also investigated the influence of selection on the posthybridization of the subspecies castaneus and musculus in Japan. Even though the genetic background of most Japanese samples closely resembles the subspecies musculus, certain genomic regions overrepresented the castaneus-like genetic components, particularly in immune-related genes. Furthermore, a large genomic block (∼2 Mbp) containing a vomeronasal/olfactory receptor gene cluster predominantly harbored castaneus-type haplotypes in the Japanese samples, highlighting the crucial role of olfaction-based recognition in shaping hybrid genomes.

Genealogical characterization of regional populations and dorsal coat color variation in the house mouse Mus musculus from Asia based on haplotype structure analysis of a gene-rich region harboring Mc1r.

We analyzed 196 haplotype sequences from a gene-rich region (250 kb) that includes Mc1r, a gene involved in coat color regulation, to gain insight into the evolution of coat color variation in subspecies of the house mouse Mus musculus. Phylogenetic networks revealed haplotype groups from the major subspecies of M. m. castaneus (CAS), M. m. domesticus (DOM), and M. m. musculus (MUS). Using haplotype sequences assigned to each of CAS and MUS through phylogenetic analysis, we proposed migration routes associated with prehistoric humans from west to east across Eurasia. Comparing nucleotide diversity among subspecies-specific haplotypes in different geographic areas showed a marked reduction during migration, particularly in MUS-derived haplotypes from Korea and Japan, suggesting intensive population bottlenecks during migration. We found that a C>T polymorphism at site 302 (c.302C>T) in the Mc1r coding region correlated with a darkening of dorsal fur color in both CAS and MUS. However, C/C homozygous mice in MUS showed marked variation in lightness, indicating the possibility of another genetic determinant that affects the lightness of dorsal fur color. Detailed sequence comparisons of haplotypes revealed that short fragments assigned to DOM were embedded in CAS-assigned fragments, indicating ancient introgression between subspecies. The estimated age of c.302C>T also supports the hypothesis that genetic interaction between subspecies occurred in ancient times. This suggests that the genome of M. musculus evolved through gene flow between subspecies over an extended period before the movement of the species in conjunction with prehistoric humans.

Whole-genome sequencing analysis of wild house mice (Mus musculus) captured in Madagascar.

In Madagascar, the house mouse (Mus musculus) is widely believed to have colonized with human activities and is now one of the most abundant rodents on the island. However, its genetic background at the genomic level remains unclear, and clarifying this would help us to infer the timing of introduction and route of migration. In this study, we determined the whole-genome sequences of five Madagascar house mice captured from an inland location in Madagascar. We examined the genetic background of samples by analyzing the mitochondrial and autosomal genomes. We confirmed that the mitochondrial genome lineages of collected samples formed a single clade placed at one of the most basal positions in the Mus musculus species. Autosomal genomic sequences revealed that these samples are most closely related to the subspecies M. m. castaneus (CAS), but also contain a genetic component of the subspecies M. m. domesticus (DOM). The signature of a strong population bottleneck 1,000-3,000 years ago was observed in both mitochondrial and autosomal genomic data. In a comparison with global samples of M. musculus, the Madagascar samples showed strong genetic affinity to many CAS samples across a wide range of Indian Ocean coastal and insular regions, with divergence time estimated as around 4,000 years ago. These findings support the proposition that the ancestors of these animals started to colonize the island with human agricultural activity and experienced a complex history during their establishment.

Insights into Mus musculus Population Structure across Eurasia Revealed by Whole-Genome Analysis.

For more than 100 years, house mice (Mus musculus) have been used as a key animal model in biomedical research. House mice are genetically diverse, yet their genetic background at the global level has not been fully understood. Previous studies have suggested that they originated in South Asia and diverged into three major subspecies, almost simultaneously, approximately 110,000-500,000 years ago; however, they have spread across the world with the migration of modern humans in prehistoric and historic times (∼10,000 years ago to the present day) and have undergone secondary contact, which has complicated the genetic landscape of wild house mice. In this study, we sequenced the whole-genome sequences of 98 wild house mice collected from Eurasia, particularly East Asia, Southeast Asia, and South Asia. Although wild house mice were found to consist of three major genetic groups corresponding to the three major subspecies, individuals representing admixtures between subspecies were more prevalent in East Asia than has been previously recognized. Furthermore, several samples exhibited an incongruent pattern of genealogies between mitochondrial and autosomal genomes. Using samples that likely retained the original genetic components of subspecies with the least admixture, we estimated the pattern and timing of divergence among the subspecies. The estimated divergence time of the three subspecies was 187,000-226,000 years ago. These results will help us to understand the genetic diversity of wild mice on a global scale, and the findings will be particularly useful in future biomedical and evolutionary studies involving laboratory mice established from such wild mice.

Phylogeographic study using nuclear genome sequences of Asip to infer the origins of ventral fur color variation in the house mouse Mus musculus.

While the house mouse (Mus musculus), widely distributed in Eurasia, is known to have substantial coat color variation between and within local populations, in both primary and secondary distribution areas, including the Japanese archipelago, the evolutionary history of the color variation is poorly understood. To address the ventral fur color variation, we quantified the lightness of museum skin specimens, and found that the southern subspecies, M. m. castaneus (CAS), has high and low lightness in dry and rainy geographic regions, respectively. The northern subspecies, M. m. musculus (MUS), has low and high levels of lightness in the high and middle latitudes of northern Eurasia, respectively. We examined sequence variation of the agouti signaling protein gene (Asip), which is known to be responsible for the ventral fur color. We performed phylogenetic analyses with 196 haplotype sequences of Asip (~180 kb) generated by phasing the whole-genome data of 98 wild mice reported previously. Network and phylogenetic tree construction revealed clustering of haplotypes representing the two subspecies, MUS and CAS. A number of subclusters with geographic affinities appeared within the subspecies clusters, in which the essential results were consistent with those reconstructed with whole mitochondrial genome data, indicating that the phased haplotype genome sequences of the nuclear genome can be a useful tool for tracing the dispersal of geographical lineages. The results of phylogeographic analysis showed that CAS mice with darker ventral fur possessed similar Asip haplotypes across the geographic distribution, suggesting that these haplotypes are major causes of the historical introduction of Asip haplotypes for darker ventral fur in mice from northern India to the peripheral areas, including the Japanese archipelago. Similarly, MUS in East Asia, which has a white abdomen, formed an Asip haplogroup with that from northern Iran, also with a white abdomen.

House mouse Mus musculus dispersal in East Eurasia inferred from 98 newly determined complete mitochondrial genome sequences.

The Eurasian house mouse Mus musculus is useful for tracing prehistorical human movement related to the spread of farming. We determined whole mitochondrial DNA (mtDNA) sequences (ca. 16,000 bp) of 98 wild-derived individuals of two subspecies, M. m. musculus (MUS) and M. m. castaneus (CAS). We revealed directional dispersals reaching as far as the Japanese Archipelago from their homelands. Our phylogenetic analysis indicated that the eastward movement of MUS was characterised by five step-wise regional extension events: (1) broad spatial expansion into eastern Europe and the western part of western China, (2) dispersal to the eastern part of western China, (3) dispersal to northern China, (4) dispersal to the Korean Peninsula and (5) colonisation and expansion in the Japanese Archipelago. These events were estimated to have occurred during the last 2000-18,000 years. The dispersal of CAS was characterised by three events: initial divergences (ca. 7000-9000 years ago) of haplogroups in northernmost China and the eastern coast of India, followed by two population expansion events that likely originated from the Yangtze River basin to broad areas of South and Southeast Asia, including Sri Lanka, Bangladesh and Indonesia (ca. 4000-6000 years ago) and to Yunnan, southern China and the Japanese Archipelago (ca. 2000-3500). This study provides a solid framework for the spatiotemporal movement of the human-associated organisms in Holocene Eastern Eurasia using whole mtDNA sequences, reliable evolutionary rates and accurate branching patterns. The information obtained here contributes to the analysis of a variety of animals and plants associated with prehistoric human migration.